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ABSTRACT 

We present an alternative method to analyse cosmic microwave background (CMB) 
maps. We base our analysis on the study of the partition function. This function is 
used to examine the CMB maps making use of the different information embedded at 
different scales and moments. 

Using the partition function in a likelihood analysis in two dimensions (Qrms-PSj"-)i 
we find the best-fitting model to the best data available at present (the COBE-HMR 4 
years data set). By means of this analysis we find a maximum in the likelihood function 
for n — 1.8tQ;g5 ^^^ Qrms-ps — IOI2.5 A^K (9^ "^ confidence level) in agreement with 
the results of other similar analyses (Smoot et al. 1994 (1 yr), Bennet et al. 1996 (4 

Also making use of the partition function we perform a multifractal analysis and 
study the possible fractal nature of the CMB sky. We find that the measure used in 
the analysis is not a fractal. 

Finally, we use the partition function for testing the statistical distribution of the 
COBE~T)MYi data set. We conclude that no evidence of non-Gaussianity can be found 
by means of this method. 

Key words: cosmic microwave background - methods: data analysis 
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1 INTRODUCTION 

In a few years the forthcoming CMB data sets from the 
missions MAP (NASA) and PLANCK (ESA) will offer us a 
much better image of the young universe than ever before. 
The CMB represents a view of the universe when it was 
about 0.002% of its present age. CMB anisotropics provide a 
link between theoretical predictions and observational data. 
Undoubtedly these data will constrain more accurately the 
fundamental cosmological parameters. In recent years sev- 
eral groups have been very active in the study of the CMB 
anisotropies. Many statistical methods have been adapted 
to the analysis of the future CMB maps, and others are be- 
ing developed. 

There are many methods which can give relatively accurate 
values for the parameters of the cosmological models. For 
example, the power spectrum is considered to be the best 
discriminator between different models (Bond et al. 1997, 
Hinshaw et al. 1996a, Wright et al. 1996, Tegmark 1996). Re- 
lated through a Legendre expansion to the power spectrum, 
the two-point correlation function is also a useful discrimina- 
tor (Cayon et al 1996, Hinshaw et al 1996b). These analyses, 
based on the power spectrum, are considered as classical but 



there are many other methods that do not make use of the 
power spectrum. The quality of the CMB maps demands for 
other statistics to supplement the power spectrum, looking 
for instance at morphological or topological characteristics 
of the data. For example the one-dimensional analysis is 
a geometrical method useful for one-dimensional scans of 
CMB data and is based on the study of regions above or 
below a certain level (Gutierrez et al. 1994) . The peak anal- 
yses, similar to the previous one but for two-dimensional 
data, deal with the number of spots above a given threshold 
(Fabbri & Torres 1995) or with other geometrical proper- 
ties like the Gaussian curvature or excentritity of the max- 
ima (Barreiro et al. 1997). Other approaches are the genus 
(Smoot et al. 1994, Torres et al. 1995), Minkowski function- 
als, which relate several geometrical aspects at the same time 
(Schmalzing & Gorski 1997, Winitzki & Kosowsky 1997), 
wavelet based techniques (Pando et al. 1998, Hobson et al. 
1998) who have shown wavelets to be very effective at detect- 
ing non-Gaussianity in the CMB, and fractality (Pompilio 
et al. 1995, De Gouveia dal Pino et al. 1995, MoUerach et 
al. 1998). 
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In this paper we present an alternative method to an- 
alyze CMB maps based on the partition function. This 
function contains useful information about the tempera- 
ture anisotropies at the different scales and moments. The 
method presented here is related to the one used by Smoot et 
al. (1994) based on moments at different smoothing angles. 
However, our method is more general and powerful because 
it works with any moment, not only with positive and inte- 
ger ones. 

The structure of the paper is as follows. In Section 2 we 
present the partition function and discuss its main charac- 
teristics. In the same section three different analyses based 
on that function are introduced: a likelihood analysis, a mul- 
tifractal analysis and a test of Gaussianity. The likelihood 
analysis uses the partition function to search for the parame- 
ters (Qrms-ps and n) that best fit a given data set. The mul- 
tifractal analysis searches for scaling laws and fractal behav- 
ior of the data. There are theoretical reasons (Sachs- Wolfe 
effect) to expect scaling in the data. Here we present the gen- 
eralized fractal dimensions and the scaling exponents and 
also comment on the possible multifractality of the CMB 
sky. In a recent paper, Ferreira et al. (1998) found evidence 
of non-Gaussianity in COBE-QMK data at 99% confidence 
level. We show that the partition function can be used to 
study this property. There is a clear relation between the 
partition function and the cumulant function, the last one 
having a specific form for a Gaussian signal. In section 3 
we apply the results of the previous sections to the COBE- 
DMR 4 years data set and we compare with other results. 
We conclude in section 4. 



2 THE PARTITION FUNCTION 

Let us start directly with the definition, 



Z{q,5) = 



Afboxc=('5) 



^^^{5Y 



(1) 



where Z{q,5) is the partition function. The quantity j-ii{5) 
is called the measure, it is a function of S which is the size or 
scale of the boxes used to cover the sample. The boxes are 
labeled by i and Nhoxcs{S) is the number of boxes (or cells) 
needed to cover the map when the grid with resolution 5 is 
used. The exponent 5 is a continuous real parameter that 
plays the role of the order of the moment of the measure. 
Let us consider a map of A'' pixels. Now the map is divided 
in boxes of size 5x5 pixels and the measure jJ.i{5) is com- 
puted in each one of the resulting boxes. Changing both, q 
and 5, one calculates the function Z{q, 5). We would like to 
emphasize that the calculation of Z{q, 5) is 0{N). 
One is free to make any choice of the measure fi{5) provided 
that several conditions are satisfied, the most restrictive be- 
ing iJ.i{5) > 0. There are no general rules to decide which 
is the best choice. For CMB maps, we use the most natural 
measure defined as follows : 



f^ii^) = Y ^^ 



pix ■ Gbox^ 



(2) 



Thus the measure in the box i is the sum of the absolute 
temperatures Tpix of the pixels inside the box in units of 
Kelvin. This is a very natural measure comparable to the 



measure used in the study of galaxies or clusters distribu- 
tion (Martinez et al. 1990, see Borgani 1995 for a review), 
where the measure is taken as the total mass (or the to- 
tal number of galaxies/clusters) contained in the box. The 
constant T* is a normalization constant. The measures are 
interpreted as probabilities and they have to be normalized, 
i.e ^. /li = 1. So T* is simply the sum of the absolute tem- 
peratures over all pixels and therefore is a constant for all 
boxes and scales. 

The temperature in the pixels is almost the same ev- 
erywhere because of the homogeneity of the signal, and one 
expects that different models will behave in a very similar 
way, making difficult the task of distinguishing them. We 
shall show how the partition function overcomes this prob- 
lem. 

Alternatively, Pompilio et al. (1995), in a multifractal anal- 
ysis of string induced CMB anisotropies (one dimensional 
scans) used as measure: 



MS) = J2 [^^ 

j=i-AIS/2 



ij + lj 



(3) 



where Aj denotes the fluctuation of the temperature in pixel 
j with respect to the mean. Here M is the total number of 
points in the data set, and (i — M5/2) and (i -I- M5/2) are 
the lower and upper edges of the ith segment with M x 5 
points, centered on the ith point of the scan. The scale 5 
runs between 1/M for the smallest segment and 1 for the 
whole segment. However, this measure is not sensitive to the 
sign of the temperature fluctuations because of the square 
in its definition. Due to this fact the full information of the 
fiuctuations is not conveniently considered. In addition, the 
generalization of this measure to 2D maps is not unique. 

Using the measure proposed in this paper, the differ- 
ences between two temperature data sets appear when high 
values of the exponent q are considered. The method is able 
to differentiate between two very close models with q rang- 
ing between [—2.5 x 10^,-1-2.5 x 10^]. This range for q is 
in agreement with the level of inhomogeneity. We are using 
absolute temperatures, that is, we have inhomogeneities of 
order 10~^ with respect to the mean value and the signal 
is almost fiat. One can consider q as a powerful microscope, 
able to enhance the smallest differences of two very simi- 
lar maps. Furthermore, 5 is a selective parameter. Choosing 
large values of q in the partition function, favors contribu- 
tions from cells with relatively high values of iJ.i{5) since 
Pi S> n'j for fii > Hj, ii q ':$> 0. Conversely, g <C favors 
the cells with relatively low values of the measure. This is 
the role played by the moments, changing q one explores the 
different parts of the measure probability distribution. The 
other parameter, 5, acts like a filter. Choosing big values of 
5 is similar to apply a large scale filter to the map. One looks 
at different scales when the parameter 5 is changed. 

To summarize, Z{q, 5) contains information at different 
scales and moments. The multi-scale information gives an 
idea of the correlations in the map, meanwhile the moments 
are sensitive to possible asymmetries in the data, as well 
as some deviations from Gaussianity. In what follows we 
show the power of the partition function to extract useful 
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information from CMB data. Three different analyses are 
used for this purpose. 

2.1 Likelihood analysis 

We shall use the partition function to encode the information 
of a given map. We compute it both for the experimental 
data and for simulated ones corresponding to different mod- 
els. In this process we are comparing the data and the models 
at several scales and using different moments. If there are 
some differences at some scale or moment, then the partition 
function should make it evident. The likelihood function will 
have a maximum for the best-fitting model to the data. For 
the CMB maps analyses, we consider models corresponding 
to different values of the spectral index n and the normal- 
ization Qrms-PS- 

The likelihood is defined in the usual way (assuming 
a Gaussian distribution for In Z{q,S)). We work with Z — 
In Z{q, 5) instead of Z{q, S) because of the large values of q 
which make imposible to compute directly Z{q, 5), 



L{Qiins-ps,n) 
where, 



(27r)"/2(deiM)i/2 '"''P^ 2 



X ) 



(4) 



X 



Y,Y.^{Z{i)) ~ Zn{i))M:,\{Z(j)) - Zn{j)), (5) 



=1 j=i 



and {Z{i)) is the average of the Z for the A'^ca realizations 
of the model at bin i. The index i defines pairs of values 
{q,5) and runs from 1 to Nq x A^^. That is, i runs from 1 to 
the total number of points Np where Z{q, 5) is defined. 
Zoii) is the value of Z for the experimental data at bin i. 



realizations: 



M,. 






{z,{i) - {z{i))){z,U) - {z{j))). 



(6) 



Zk{i) denotes the value of Z at bin i for the k realization. 
We tried different number of realizations A'^roa but the results 
appear to be stable for A'^roa > 2000 per value of Qrms-ps 
and n. 

We have two possibilities to perform a best fit to the data. 
The first one is to minimize x^ a-nd take the values of 
Qrms-PS and n at the minimum of the x^ surface as the 
best-fitting values. The second possibility is to work with 
the likelihood L looking for the maximum. We tested the 
two possibilities using simulated CMB maps derived from 
a given pair of parameters (Qrms-PS, n) and then using 
these maps as the data maps. Due to cosmic variance we 
obtain a set of maxima in the likelihood and of minima in 
the x^- The conclusion is that the likelihood is somewhat 
better than the x^ as expected. For instance with 2000 in- 
put realizations with Qrms-PS = 14 fiK. and n = 1.3 the 
distribution of maximun likelihood values finds a maximum 
at Qrms-PS = I3I4 25 mK and n = 1.2'^°q\^ while the x^ ren- 
ders a minimum in Qrms-PS = 17^4 /.iK and n = l.Olg^g. 
The errors are marginalised at the 68% confidence level and 
are similar to those obtained with the standard methods 
based on the power spectrum (see for instance Wright et al. 
1996). 



2.2 Multifractal analysis 

The notion of multifractal measure was first introduced by 
Mandelbrot (Mandelbrot 1974) in order to study different 
aspects of the intermittency of turbulence (see also Sreeni- 
vasan and Menevau 1988). The multifractal formalism was 
further developed by many other authors and today it is a 
standard tool applied in almost all fields of science: molecu- 
lar physics, biology, geology, astronomy, etc . In the context 
of the description of the large scale structure of the Universe 
it was first introduced by Jones et al. (1988). 

Some authors (Pietronero & Sylos Labini 1997) sug- 
gest that the distribution of matter in the universe is fractal 
with dimensionality D2 — 2. They defend that the scal- 
ing remains up to the larger scales probed by the present 
day available redshift catalogues. Many other authors, how- 
ever, have found enough evidence of homogeneity at large 
scales (Davis 1997, Guzzo 1997, Scaramella et al. 1998) in 
the analysis of the same data sets. One of the basic tenets of 
the standard cosmology is that at very large scales the dis- 
tribution of matter is homogeneous. The homogeneity and 
isotropy of the CMB support this overwhelming evidence, 
indicating that there exists a continuous transition between 
scale invariant clustering at small scales and homogeneity at 
large scales (Martinez et al. 1998; Wu, Lahav & Rees 1998). 

At large angular scales, the CMB anisotropies AT/T 
generated from a scale-free density perturbation power spec- 
trum in a flat Q. = 1 universe can be described by a frac- 
tional Brownian fractal (as shown in MoUerach et al. 1998). 
In particular, both inflationary and defect models predict an 
approximately scale invariant Harrison-ZeVdovich spectrum 
on large angular scales showing the scaling predicted by the 
Sachs-Wolfe effect. At small angular scales (0.2° < 9 < 1°) 
the predictions of inflation and topological defects models 
are different (Durrer et al 1997) allowing to differentiate 
them. It is then interesting to study the possible fractality 
of the CMB anisotropies, since the seeds or fluctuations that 
are supposed to be the precursors of the largest structures 
observed today, are yet unperturbed by evolutionary phe- 
nomena. Several works follow this kind of analysis. In the 
paper by De Gouveia dal Pino et al. (1995), the authors 
based their analysis in the study of the perimeter-area rela- 
tion of the isocontours of temperature at a given threshold. 
They used the COBE-QMR 1 year data set and only the 53 
GHz channel. They found evidence for a fractal structure in 
the COBE-QMK data with dimension D — 1.43 suggesting 
that the CMB could not be homogeneous. Apart from the 
fact that these data have a low signal to noise ratio, this 
does not necessarily mean that the CMB is not homoge- 
neous. This dimension corresponds to the temperature iso- 
contours, and not to the temperature itself. Other works use 
multifractal analysis with CMB. Pompilio et al. (1995) apply 
the multifractal analysis to simulated string-induced CMB 
scans searching for the non-Gaussian behavior induced by 
cosmic strings. More recently, MoUerach et al. (1998) have 
applied a fractal analysis in order to study the roughness of 
the last scattering surface and used this technique to search 
for the model that best fit the COBE-DMR 4yr data. These 
authors show the capabilities of this method for the analysis 
of future data, in particular for those experiments with high 
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signal to noise ratio. 

In this section we will use the partition function to study 
the possible multifractality of the CMB sky, using as mea- 
sure the absolute temperature (see eq. 2). The multifrac- 
tal analysis has been presented in several versions but the 
most popular is due to Frisch and Parisi (1985), Jensen et 
al. (1985) and Halsey et al. (1996), where the spectrum of 
singularities /(a) was introduced. We will give here a brief 
description of the multifractal approach. A presentation of 
the method can be found in Feder 1988, Schuster 1989, Vic- 
sek 1989 and more formally in Falconer 1990. 
The multifractal formalism has as starting point the parti- 
tion function. The generalized or Renyi dimensions are de- 
fined by the asymptotic behavior (as the scale 5 tends to 
zero) of the ratio between \nZ{q,5) and In 5, 



i^(,) = lim^^"^(f). 

It is easy to see that for q - 
or capacity dimension, 

lnA'"boxcs(5) 



D(0) = lim 



s^o ln(l/(5) 



(7) 
wc obtain the box-counting 



(8) 



For q — 1, D{1) is the information dimension, which is ob- 
tained from Eq. (7) by applying L'Hopital's rule. For q — 2, 
D{2) is the correlation dimension (see Schuster 1998 for 
other alternative definitions and the relation between them). 
A simple fractal or monofractal is defined by a constant 
D{q). Dependence of _D on g defines a multifractal. In most 
of the practical applications of the multifractal analysis, the 
limit in Eq. (7) cannot be calculated, either because we do 
not have information for small distances (as it happens in 
this case) or because below a minimum physical length no 
scaling can exist at all (for example the size of a galaxy in the 
multifractal nature of the galaxy distribution). This problem 
is usually overcome by finding a scaling range [5i, 82] where 
a power-law can be fitted to the behavior of the partition 
function 



Z{q, S) oc 5"<«' for Si<5<S2 



(9) 



The scaling exponents T{q) are related with the generalized 
dimensions by 



r(g) = (q-l)D(g). 



(10) 



Other quantity, commonly used in the characterization of 
multifractals, is the so-called /(a) spectrum. If for a given 
box (labeled by j) the measure scales as 



w('5)~^"^ 



(11) 



then, the exponent a, which depends in principle on the 
position is known as crowding index or Holder exponent. If 
all the points have the same scaling, then all the exponents 
a will be the same and this corresponds to a monofractal. 
Otherwise, if we have boxes with different scaling, what we 
have is a mixture of monofractals. This set is known as a 
multifractal (each monofractal formed by the points with 
the same scaling and therefore with the same exponent a). 
The exponent a is used to label the boxes covering the set 
supporting a measure, thereby allowing a separate counting 
for each value of a. In a multifractal set a can take different 
values within a certain range, corresponding to the different 



strength of the measure (Halsey et al. 1996). The subset 
formed by the boxes with the same a will be denoted 5*0,. 
This subset has Nc,{5) elements (boxes) and in general, for 
a multifractal set, this number varies with the scale S as 



N^{5) r^S-f-^K 



(12) 



Comparing this expression with the definition of the box- 
counting dimension, Eq. (8), the quantity f{a) can be inter- 
preted as the fractal dimension of the subset Sa- However, 
this physical meaning of the function /(a) is not always true 
(Grassberger, Badu & Politi 1988; Falconer 1990). 

It can be shown (Halsey et al. 1996; Martinez et al. 
1990) that the quantities q and r(q) can be related through 
a Legendre transformation with a and /(a). These relations 



a{q) = 



dr{q) 



f[a) = qa{q) -riq). 



(13) 
(14) 



To illustrate this section we use the well known multi- 
plicative multifractal cascade (Meakin 1987; Martinez et al. 
1990). The construction of this multifractal is as follows: A 
square is divided into four equal square pieces and a prob- 
ability Pi, (i = 1, . . . , 4), such that '^ _-.Pi = 1, is assigned 
to each one. Each piece is again subdivided in four small 
squares, allocating again a value pi randomly permuted to 
each one. The measure assigned to each one of the new sub- 
squares is the product of this value of pi and the corre- 
sponding value of its parent square. The subdivision process 
is continued recursively. In Fig. 1 we show a realization of 
this multifractal on a grid of 256 x 256 pixels for the values 
of the probabilities pi — 0.18, P2 ~ 0.23, ps = 0.28 and 
P4 = 0.31. We can easily calculate the theoretical values of 
the multifractal functions D{q) and /(a) for this illustra- 
tive example (Martinez et al. 1990) . With the multiplicative 
multifractal we tested the power of the method to recover 
the true dimensions. In Fig. 2 we show the generalized di- 
mensions D{q) and the corresponding spectrum of fractal 
dimensions f{a). These curves match perfectly the theoret- 
ically expected ones. Note that a single monofractal should 
render a straight line for D{q) and a single point for /(a). 



2.3 Testing Gaussianity 

A Gaussian distribution of CMB temperature fluctuations 
is a generic prediction of inflation. Forthcoming high- 
resolution maps of the CMB will allow detailed tests of 
Gaussianity down to small angular scales, providing a crucial 
test of inflation. Most of the works that analyse CMB maps 
assume Gaussian initial fluctuations. Kogut et al. (1996) 
find that the genus, three-point correlation function, and 
two-point correlation function of temperature maxima and 
minima are all in good agreement with the hypothesis that 
the CMB anisotropy on angular scales larger than 7° rep- 
resents a random-phase Gaussian field. Other alternative 
methods are proposed, like the angular- Fourier transform 
(Lewin et al. 1998), Minkowsky functionals (Schmalzing & 
Gorsky 1998), correlation of excursion sets (Barreiro et al. 
1998), and the bispectrum (Heavens 1998). In an analysis 
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of the 4 years COBE-DMR data based on the bispectrum 
Ferreira et al. (1998) have found that Gaussianity is ruled 
out at a confidence level in excess of 99% near the multipole 
of order I = 16. 

In this section we will test the Gaussianity of the CMB 
data using an alternative method. The idea is to use the 
relation between the partition function and the generating 
function, the last one defined as, 



G4t) = (e'l- 



(15) 



If we know that x is Gaussian distributed then, solving the 
integral corresponding to the mean value of the previous 
definition, results : 



{t) = e 



t{x) + - 



(16) 



It follows from the last expression that the cumulant func- 
tion F is: 



F{t)=lnG{t)=t{x) + ^. 
Finally, the function 

2 2 

H{t) ^ F{t) ~ t{x) ~ ^-^ 



(17) 



(18) 



should be zero for all t for a Gaussian field. 

Let us consider the definition of Z{q,S). If the measure is 

defined as 



M*(5) = e^'W, 



(19) 



with iJ.i{S) the same as in section 2. Then the partition func- 
tion is 



sW 



Z{q,5) 



E 
1=1 



o9Mi(*) 



or equivalently, 

Z{q,5) = iVboxcs(<5){e'^) = iVboxcs(5)Gp(g). 



(20) 



(21) 



This relation between Z{q,5) and G^{q) allows us to con- 
struct the function H{q) which, for a Gaussian measure fi, 
should be zero for all q at each scale S. This is a simple 
way to find non-Gaussian signals. The function H{q) rep- 
resents the contribution of all the moments larger than 2. 
This contribution should be zero only for a Gaussian field. 
A plot of this function indicates directly the deviations from 
Gaussianity. 



3 RESULTS: APPLICATION TO COBE-DMR. 
DATA 

As a practical use of the methods presented we will apply 
them to the 4 years COBE-DMR data. 



3.1 Description of the data 

We use the COBE-DMR 4 years 53 + 90 GHz maps combi- 
nation, which is the choice with the largest signal to noise 
ratio (Bennet et al. 1996). These data are in the Quad- Cube 
pixelization with a pixel size of ~ 2.6° and the resulting 
number of pixels is 6144. The data in each pixel represents 



AT in mK units. The dipole has already been subtracted. 
Assigned to each pixel there is an additional information, 
the number of times that this part of the sky was explored 
by the antenna. This information is relevant for the estima- 
tion of the instrumental noise. 

Part of the data is contaminated by the galactic emission. 
There is a strip between ~ ±20° (in galactic coordinates) 
in which the galactic emission dominates the CMB signal. 
This strip should not be included in the analysis in order 
to avoid spurious signals. In addition to this strip there are 
two patches in the sky (one near Orion and the other one 
in Ophiucus) that show a strong galactic emission at mm 
wavelengths (Cayon et al. 1995), and should therefore also 
be removed from the analysis. When this mask is applied, 
the number of surviving pixels reduces to 3881 from the 
original 6144. 



3.2 Likelihood analysis 

In order to determine which are the values of the quadrupole 
normalization Qrms-ps and the spectral index n that best 
fit the COBE data, we perform Monte Carlo simulations of 
the CMB maps for a scale-free model with a power spec- 
trum given by P{k) oc k", which has variance in the aim 
multipoles given by (Bond and Efstathiou 1987): 



C, 



5^ m 



r[/ + (n-l)/2]r[(9-ri)/2] 
'r[l + {5-n)/2]r[{3 + n)/2Y 



(22) 



We consider difi'erent values for Qrms-ps and the n rang- 
ing from Qrms-PS = 4/iK to Qrms-PS = 35/iK and from 
n = 0.3 to n = 2.3. We add instrumental noise based on 
the number of data collected by COBE-DMR at each pixel. 
Furthermore, there is another effect that must be taken into 
account, the cosmic variance. To treat conveniently this ef- 
fect we perform a large number of simulations (> 2000) for 
each pair of values (Qrms-ps , n) and then we compare the 
average Z = \n Z{q, S) values of these simulations with the 
Z corresponding to the COBE-DMR data (the used values 
for q and 5 were q = -120000,-40000,72000,152000 and 
5 = 3,4, 8, 16 pixels). The choice for q and 5 values is based 
on the test described in the last part of section 2.1. The size 
of the Z{q, 5) grid, A^, x Ng is not critical and what is now 
relevant is the q values considered. In particular, high order 
moments (i.e large q) are very sensitive to the tail of the dis- 
tribution and therefore the results obtained with those high 
values on the parameter estimates are not stable. The com- 
bination of q and 5 values, was one of the combinations for 
wich the recovered parameters Qrms and n were closer to the 
input parameters and with smaller error bars. As mentioned 
in section 2, q should take values of order lO'' in order to 
distinguish between models with temperature fluctuations of 
order 10""". The values of q where chosen to be asymmetric 
in an attempt to consider possible asymmetries that could 
exist between the negative and positive temperature fluctu- 
ations. The range for 5 runs from 2 pixels (approximately 
the antenna size) to 24 pixels which is the largest box size 
required to have at least 8 boxes. Using a maximum likeli- 
hood method one can determine which are the best-fitting 
parameter values of the simulations (signal -I- noise) to the 
COBE-DMR data. 
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In Fig. 3 we show a contour plot of the hkeUhood 
obtained for the COBE-DMK data. The maximum is at 



Qr 



10 



+3 



fjK. and n = 1. 



(95% marginalised 



errors) and the contour level at 68% is compatible with 
the assumed standard value Qrms-ps = 18 ± 3 /iK for the 
Einstein-de Sitter model with a scale invariant primordial 
spectrum of density perturbations, n = 1. The various anal- 
ysis of the 4 years COBE data combined give as the best- 
fitting parameters Qrms-PS ~ 15.3l2'g mK and n — 1.2±0.3. 
The result presented here predicts larger values of n and 
smaller values of Qrms-PS than the result indicated above 
(although always inside the anticorrelation law for the two 
parameters). This result is in agreement with the one found 
by Smoot et al. (1994), using a similar approach. Smoot et 
al. (1994) found for the best fit, Qrms-PS = 13.2 ± 2.5 ^K 
and n = l.Tlgg. A possible explanation for the discrepancy 
between our results and those obtained with the standard 
methods could be a bias present in the likelihood estimator. 
In the tests of our algorithm we found a systematic bias in 
the marginalized likelihood functions both for Qrms and n 
with tipical values of 5n ~ -1-0.2 and 5Qrms ~ —2 which 
could explain part of our discrepancy. The reason for this 
bias can be the difference between the assumed Gaussian 
form for the likelihood of the partition function in eq. (4) 
and the real non-Gaussian distribution. The probability dis- 
tribution of the Z at each (g, 5) obtained from simulations 
is similar to a Gaussian probability distribution but with a 
longer tail for high values. We also think that maybe the 
noise can contribute to that bias. The high order moments 
(large q) of the partition function are very sensitive to the 
tails of the distribution of the temperature fluctuations. A 
low signal to noise ratio (as is the case for the COBE-DMR 
data) could raise the parameter n that best fit the COBE- 
DMR data. We did some tests in this direction and appar- 
ently the noise can increase the value of n (and consequently 
can produce a lower value of Qrms)- 



3.3 Multifractal analysis 

We apply the formalism of section 2.2 to the simulations and 
to the COBE-'DMR data. In Fig. 4 we plot D(q) and /(a) 
for the COBE-DMR and for one model (Qrma-ps = 15mK, 
n = 1.2) inside the QrmB-PS-n degeneration with its error 
bars. The D[q) curve has been obtained by fitting a power- 
law to the partition function in the range of scales 2 < 5 < 24 
pixels, following Eqs. (9) and (10) . Note that the value of 
D(0) is not 2 as it would be expected for a continuous bidi- 
mensional surface. The mask slightly lowers this value. 
By means of a Legendre transform (Eqs. (13) and (14)) we 
have obtained the corresponding /(a) curve. A narrow /(a) 
curve means a very homogeneous data. If the measure as- 
sociated to the data is multifractal in nature, these curves 
should be the same for all the scale ranges. We have found 
that this is not the case for the COBE-DMR data. The 
multifractal curves corresponding to difi'erent scale ranges 
do not match each other. CMB simulations without noise 
show the same behavior. The reason for that lies in the fact 
that a scaling like that in Eq. (9) is not present. This can 
be illustrated by looking at the behavior of the local slopes 
of \nZ{q,5) vs. In (5. In Fig. 5 we show the change in the 
reduced slopes {r[q, S)/(q — 1)) as a function of the scale for 
a fixed value of the parameter q for the COBE-DMK data. 



For this plot the analysis was performed only in the top and 
bottom faces of the Quad-Cube which are not affected by 
the mask. For a multifractal measure these curves should 
be horizontal straight lines. As we can appreciate in the left 
panel, this is not the case for the COBE-DMK data. The 
result for a simulation without noise is shown in the right 
panel. In both cases, we do not see a neat plateau for large 
absolute values of q. However it is not clear whether the 
fiuctuations of the local reduced slopes are just due to nu- 
merical noise related to the resolution of the maps (i.e. the 
limited number of pixels) or, on the contrary, these fluc- 
tuations are intrinsic to the measure and, therefore, prove 
that the measure is not a multifractal. Although our result 
neither support nor contradict this interpretation, it seems 
more natural to expect fractal behaviour in the case that 
one is using the absolute value of the relative temperature 
fiuctuations AT/T as the measure (MoUerach et al. 98). As 
shown in that paper, AT/T fiuctuations generated by the 
Sachs-Wolfe effect behave like a fractional brownian fractal. 



3.4 Gaussianity 

To test whether the COBE-DMK data are Gaussian dis- 
tributed, we compare the H{q) curve for COBEJ-DMK with 
those curves arising from the best-fitting CMB Gaussian 
models obtained in section 3.2. In Fig. 6 we show the plots of 
H{q) for different grid scales. For each realization, th mea- 
sure is rescaled in order to have dispersion equal to one. 
This allows to have a small and equal range of q values for 
all scales. We would like to point out the deviation of the 
mean value from zero when q moves away from zero. This is 
due to the fact that we have a finite number of pixels (i.e. a 
cosmic variance effect). The predicted behavior of equation 
(15) is only true when we compute the mean over infinite 
values (or equivalently, solve the integral between — oo and 
-l-oo). Otherwise H{q) is not zero at large values (positive 
and negative) of q. Fig. 7 shows the likelihood distribution 
of the 1000 Gaussian realizations (with noise) and the dot- 
ted line corresponds to the COBE value. It is clear that 
the COBE-DMK is perfectly compatible with the Gaussian 
hypothesis. 



4 DISCUSSION AND CONCLUSIONS 

We have shown in this work the power of the partition func- 
tion to describe CMB maps taking into account the infor- 
mation given at different scales and by different moments. 
We have also shown the flexibility of such a function to 
be used in various analyses: standard likelihood, multifrac- 
tal and Gaussian. We applied these analyses to the 4 years 
COBE-DMK data. 

Based on the likelihood function we flnd the best-fitting 
parameters Qrms-PS ~ IOI2 5 mK and n = l.S^Qgg. It is re- 
markable the agreement between our work and the one by 
Smoot et al. (1994). 

The COBE-DMK data (and the simulations of scale invari- 
ant power spectrum) do not show a fractal behavior, regard- 
ing the absolute temperature map. On the other hand, re- 
cent galaxy surveys covering large scales (> 100 Mpc) do not 
show either a fractal behaviour (Wu, Lahav & Rees 1998). 
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Both results allow to conclude that neither the mass distri- 
bution (assuming a linear bias) nor the intensity of the CMB 
show a fractal behaviour on large scales. 

The partition function analysis performed shows no ev- 
idence for non-Gaussianity in the COBE-DMR data. This 
is in agreement with all the previous analyses of the COBE- 
DMR data except the one by Ferreira et al. (1998). Simula- 
tions done at higher resolution have shown the power of this 
method to discriminate between Gaussian and non-Gaussian 
signals. That analysis will be presented in a future paper. 

Finally, we would like to remark that the likelihood 
analysis based on the partition function is computation- 
ally intensive. Actually a non optimized code applied to the 
COBE-DMR data takes a few days (CPU time) to run in 
an Alpha server 2100 5/250. Moreover, the computation of 
the partition function increases with the number of pixels 
N as 0{N). This rate should be compared with the most 
widely applied method used to compress data and to esti- 
mate cosmological parameters, the power spectrum of the 
fluctuations. The direct computation of the power spectrum 
goes like 0(A''logA'') (this behaviour is due to the FFT). 
Standard brute- force approaches used to estimate the power 
spectrum go like 0{N^). The reason for this 0{N^) rate is 
the matrix inversion and determinant calculation whose di- 
mension grows as the number of pixels. On the contrary, 
in the partition function likelihood analysis the number of 
bins (or matrix dimension) of the likelihood is A''^ x Ng, be- 
ing this number usually well below one thousand (even for 
high resolution maps). The number of moments q is an arbi- 
trary parameter independent of N and the number of scales 
5 increases as < 0(A^"). The process of inverting the cor- 
relation matrix is clearly reduced in the case of the partition 
function. This point makes the method useful for forthcom- 
ing large data-sets. One can therefore consider the partition 
function as an alternative way to compress large data sets. 
Furthermore, for the general situation of non-Gaussian data 
sets, the partition function is clearly preferable to the power 
spectrum since the former contains information on several 
moments of the data. 
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Figure 2. D(q) (left) and f(a) (right) for the multiplicative mul- 
tifractal. The dotted line on the right panel represents the theo- 
retical /(a) curve for a larger range of q 



Figure 3. Contour confidence levels (68%, 95%, 99%) for the 
likelihood distribution. The maximum is at Qrms = 10 /xK and 
n = 1.8. 



Figure 4. D{q) (bottom) and /(o) (top) curves for COBE-'DMK 
(dotted line) and for one model with Qrms = 15 /iK and n = 1.2. 



Figure 5. Reduced slopes computed between two consecutive 
scales for the COiJi?-DMR data (left) and for simulated data 
(without noise, right). Each curve corresponds to a fixed value of 
q betwen — le6 and le6. 



Figure 6. H{q) curves for eight different scales. From left to right 
and from top to bottom. 2x2 pixels to 24 X 24 pixels. The x-axis 
represents the values of q. 
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Figure 7. Likelihood distribution for 1000 simulations for the 
best-fitting model in section 3.2. The value corresponding to the 
COBE-DMR data is shown as a dotted line and is in clear agree- 
ment with the Gaussian hypothesis. 



© 0000 RAS, MNRAS 000, 



2.8 



2.6 



2.4 - 



2.2 






2.0 



Generalized fractal dimensions 




50 



100 

q 



150 



200 



Spectrum of fractal dimensions 



^ 1.0 - 




40 



30 






20 - 



Ikellhood 



o 



10 




_J \ \ \ I \ \ I \ \ I 

1 2 

spectral index n 



cd 
Si 




0.5 



-0.5 



1.93394 



1.93395 



alpha 



1.93396 



1.93397 



1.93397 



1.93396 



Q 



1.93395 



1.93394 




-4x10 



4x10° 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



2.00004 



2.00002 



CD 

o 

CO 



1 .99998 



1 .99996 




T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



2.00004 



Model (without noise) 



2.00002 



q = - 1 000000 



CD 

o 

CO 



1 .99998 




q = + 1 000000 



1 .99996 



J I I I I I I I I I I I I I I I I I I I I I I I I I I L 



3 4 
Scale 




T 

7 


M M M M 


1 r 


' 1 


1 hill 


^ Iji 


II 1 1 r 

CO 


' 


4 :[ 


— 


- 


1 1 1 1 1 1 1 1 1 t 



-5 5 



-5 5 







-8 








-2 -r 



-8 




-5 5 



-5 5 













-5 5 



-5 5 





-5 5 



-5 5 



T 1 1 r 



100 



80 



60 



40 



20 



T 1 1 r 



~i 1 n 



m 



~i n_ 



_n n_ 




J I I L_ 



J I I L 



J L 



Likelihood 



